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Colour octet contribution to exclusive P-wave charmonium decay into 
nucleon-antinucleon 

S.M.H. Wong * a 

a Nuclear and Particle Physics Section, University of Athens, Panepistimiopolis, GR-15771 Athens, 
Greece 

Institute of Accelerating Systems and Applications (IASA), P.O. Box 17214, GR-10024 Athens, Greece 
and 

Fachbereich Physik, Universitat Wuppertal, D-42097 Wuppertal, Germany 

We show that although there is no infrared divergence in exclusive P-wave quarkonium decays, the colour 
octet contribution is no less important than in the inclusive decay. Results from more proper calculations with 
only colour singlet contribution are not sufficient to account for the measured partial decay widths and so the 
next higher Fock state must be included. Using the \j decay into NN as an example, a scheme is devised 
to yield estimates of the decay widths using both contributions, the results are comparable with experimental 
measurements. IASA 98-2, UA/NPPS-98-12 
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1. Introduction 

In the late 60's and early 70's, Barbieri et 
al |l],H considered the inclusive decay of P- 
wave heavy quark-antiquark system into hadrons. 
Such decays proceed via the annihilation of the 
heavy fermions into massless gluons. For P-wave 
quarkonium, annihilation into two gluons is al- 
lowed and this is the leading contribution to the 
inclusive hadronic decay width, at least for the 
total angular momentum J = and J = 2 sys- 
tem. For J — 1 P-wave quarkonium, the decay 
into two massless spin-1 gluons is forbidden. The 
leading decay process in this case is therefore one 
gluon and a light quark-antiquark pair. Whereas 
the kinematics of the leading decay process of the 
even- J systems into two gluons is very simple and 
completely well defined, that of the J = 1 system 
is not. In the rest frame of the quarkonium, the 
decay into two gluons would proceed with them 
moving apart from each other with equal and op- 
posite momentum. The three-body leading de- 
cay of the J — 1 quarkonium is, on the other 
hand, kinematically less restrictive so the situa- 
tion where the quark and antiquark move apart 
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with nearly equal and opposite momentum with 
an accompanying very soft gluon is allowed. The 
corresponding probability on the tree graph level 
from such a process with a very soft gluon is easily 
shown to contain infrared divergence when all the 
participants are on their respective mass shells. 
In view of the infrared divergence and the fact 
that the heavy quark and antiquark are really in 
a bounded system, one should perhaps keep them 
off shell by an amount corresponding to the bind- 
ing energy e of the system. Thus using e as an 
infrared cutoff, one gets for the J = 1 quarkonium 
decay width 

dq 

q 

This is the solution proposed in In fact, 

this infrared divergence also appears in the de- 
cay widths of the even- J P-wave quarkonia but 
at the next-to-leading order. 

With the more recent advances in quarkonium 
physics we now know that this infrared di- 

vergence should not be present. The physical pic- 
ture of the quarkonium in the early 70's was a 
heavy quark bounded with an equally heavy an- 
tiquark. The missing component in this picture 
is the higher Fock states of the quarkonium. The 
contribution to the decay into gqq from the so- 
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called colour octet, the next higher Fock state 
of the P-wave quarkonium, becomes degenerate 
with that from the lowest valence state when the 
gluon approaches the infrared. When both are 
included, there is no infrared divergence. This is 
the situation of the inclusive P-wave quarkonium 
decay. Colour octet is needed and is introduced 
to cancel infrared divergence. 

In exclusive decay, however, the heavy quark 
system decays from a bounded system into other 
bounded systems, the same infrared divergence in 
the inclusive decay does not appear. So there is 
no need for colour octet and the partial widths 
of the various exclusive decay modes have been 
calculated from the valence colour singlet compo- 
nent alone, see J||| for example. We will show 
that, on the contrary, the colour octet is very im- 
portant for P-wave quarkonium. 

2. Angular momentum suppression of P- 
wave wavefunctions 

The decay of a quarkonium of mass M through 
annihilation is a short distance process with the 
annihilation length given by L ~ 1/M. For a 
heavy quarkonium, 1/M is very small so we only 
need the quarkonium wavcfunction near the ori- 
gin ip(L ~ 0). This is true only for a S-wave 
quarkonium, for a P-wave, the wavcfunction at 
the origin is zero so one expands the wavcfunction 
around the orgin and uses instead the quantity 
Lip 1 (L ~ 0) for the infrared confinement physics. 
In momentum space, this means 



S-wave: 
P-wave: 
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Assuming that there is no significant difference 
between the S- and P-wave Fourier components 
i/j(k), which seems reasonable, it follows that the 
P-wave wavefunction is weighed down by 1/M in 
comparison with that of the S-wave. Therefore 
P-wave is relatively suppressed on the level of the 
wavefunction. We will see that this is very im- 
portant in the next section. 



3. Comparing colour singlet with octet 

In order to find out whether the colour octet 
component in exclusive quarkonium decay is neg- 
ligible or not when compared with the colour sin- 
glet, we now perform power counting in the only 
large scale of the process, namely, M. We must 
choose an explicit process and a scheme for the 
calculation of the decay probability amplitude in 
order to be able to do this. We use \J decay 
into NN as the example. For the purpose of 
power counting, it is sufficient to choose the sim- 
pler standard hard scattering scheme of Brodsky 
and Lepage Jjj. In this scheme, the probability 
amplitude is given by a convolution of distribu- 
tion amplitudes <f>{x) and hard perturbative part 

and this is related to the partial decay width by 
r ~ \A4\ 2 /M. So M. must have mass dimension 
one. On the right hand side of Eq. (Q), the only 
dimensional quantities are the decay constants /'s 
and a certain power of 1/M hidden in Th- They 
must make up the right dimension to match the 
left hand side. The decay constants have the fol- 
lowing mass dimension, fjsr, fpj and the octet 
constant are from 3-particle wavefunctions and so 
must have mass dimension two. The singlet con- 
stant f®j is from a 2-particle wavefunction which 
would mean dimension one but being a P-wave 
increased this to two. Now extracting sufficient 
power of 1/M from Th to make up the right di- 
mension for M. gives 
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We see that both the colour singlet and octet are 
weighed by the inverse fifth power of the large 
scale M. Therefore the colour octet is not sup- 
pressed when compared to the singlet contribu- 
tion. 

Now we are ready to return to the 1/M sup- 
pression of P-wave wavefunction discussed in the 
previous section. Only because of this suppres- 
sion is the singlet amplitude scale like 1/M 5 . If 
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this were not present as in the case of the S-wave 
J/tp, the singlet amplitude would scale like 1/M 4 
instead. In that case, it would be legitimate to 
neglect the colour octet contribution. 

4. Calculation — singlet contribution 

To obtain the colour singlet contribution to the 
decay width, instead of the original hard scat- 
tering scheme, we employ the modified version 
of Botts, Li and Sterman [SB . In this modified 
scheme, the probability amplitude is now a con- 
volution of wavefunctions, hard perturbative part 
Tjj and the exponential Sudakov suppression fac- 
tor 

M ~ tp XJ (x,k±) ® ipN(x,kx) 8>i>N(x,k±) 

®f H {x, k±,a s (x)) <g> exp{-S{x, fcj_)} . (7) 

The essential differences between this and the 
standard scheme are that internal transverse mo- 
menta are kept everywhere, the strong coupling 
a s is part of the convolution instead of a constant 
and Sudakov suppression is included. Actually 
even in the standard scheme, a s need not be a 
constant and can be made part of the convolution, 
but in that case, some adhoc arrangement such as 
freezing the value of the coupling at small momen- 
tum transfer will have to be made. The modified 
scheme has no such problem as the more natural 
Sudakov factor takes care of potential divergence 
of the coupling. This is achieved by using the ad- 
jacent largest virtuality in the neighbourhood of 
each coupling for its argument JT^JT^] . It can be 
that of the neighbouring gluon propagator, quark 
propagator or the inverse internal transverse sep- 
aration of quark pairs in the outgoing nucleon- 
antinucleon. In this way, the scale of the decay 
process is determined dynamically. 

Our objection to some of the calculations based 
on using the standard scheme is that if a fixed 
coupling is not used then freezing of the coupling 
or one of the other methods to deal with this is 
necessary. These methods have not a sound phys- 
ical basis. If a constant coupling is used, then tak- 
ing our current example of xj decay into nucleon- 
antinucleon the probability amplitude goes like 
M. ~ o? s and the decay width is therefore r ~ a®. 
A choice of a s = 0.3 instead of a s — 0.5 will 



therefore make a factor of 20 difference in the 
width so provided the singlet countribution is of 
reasonable size, one can always argue for a value 
of the coupling so that the singlet contribution 
is sufficient to account for the measured partial 
width. This in our opinion is quite arbitrary and 
unsatisfactory 

Using the modified scheme, we found that the 
colour singlet component makes up of no more 
than 6 and 12 % of the measured decay widths 
for x\ an d Xi respectively. So colour singlet alone 
is not sufficent in agreement with our theoretical 
expectation. 



5. Calculation — including octet contribu- 
tion 

To complete the calculation, one must include 
the colour octet contribution. It turned out that 
to use the modified scheme to include also the 
octet contribution was rather daunting because 
the internal transverse momentum must be kept 
everywhere and with an extra constituent gluon, 
the dimension of the numerical integrations be- 
came quite high. One could go back to the stan- 
dard scheme but then had to face the objections 
raised in the previous section. Since the prob- 
lem with using the standard scheme is mainly in 
the coupling, one can circumvent this by using a 
model coupling that is free of the Landau pole 
and thus can be included in the convolution of 
probability distribution amplitudes. With such a 
model, the XJ decay probability amplitude into 
NN becomes 

M ~ fxAxj{ x ) ® fN<f>N(x) ® ffr<pfi(x) 

^(^^'(i)). (8) 

There exist various models that fit our require- 
ment. We choose the one constructed by Shirkov 
and Solovtsov |lo) which is less complicated and 
the Landau pole is removed by simple subtrac- 
tion. To one-loop, it is 

, , , . 47T ( 1 A 2 i 
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Table 1 

Prelimary estimates of the colour singlet plus 
octet contributions to \J decay into NN in our 
scheme. 



J 


T(xj NN) [eV] 


PDG [eV] 


1 


89.0 


75.68 


2 


180.0 


200.00 



Although this model is comparatively simple, it 
agrees with known estimates of the quantity 

MQ) = ^ [ dka s (k) (10) 
W Jo 

from jet physics. These estimates are done at a 
value of Q ~ 2 GeV 

Without the fc^ dependence, the calculation 
becomes much simpler and with our so con- 
structed semi-modified hard scattering scheme, 
the dynamical setting of scale by the process it- 
self is preserved. Details of our calculation par- 
allel those of P|| . Our preliminary results of 
the total colour singlet plus octet contribution to- 
gether with the experimental measurements fl~5| 
are shown in Table 1 jl6|. As can be seen, the 
agreement with the measured values are quite rea- 
sonable. This confirm our theoretical argument in 
Sec. [| I n fact, this argument can be extended to 
even higher wave quarkonia and in such cases, not 
only the next higher Fock state must be included 
but also the next-next higher states. 

In view of the above considerations and results, 
the constituent gluon in P-wave quarkonia is an 
important part of the heavy quark hadrons. The 
description that they are bounded heavy quark- 
antiquark systems becomes inaccurate. If the 



"quark" in quarkonium means heavy quark- 
antiquark object, then the three Xj' s would each 
be less of a quarkonium than a J/ij) and a D-wave 
would be even less so. Therefore the nomencla- 
ture is missleading in this sense. 
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